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ABSTRACT 

I show that the cluster Lx — o~ relation should be sensitive to cosmologies with a scale- free power spec- 
trum of initial density fluctuations, P(k) oc k n . I derive the dependence, and argue that a conservative 
interpretation of current observations implies n < —2.0 and n < —1.1 at the one-sided 90% and 99% 
confidence levels, respectively. This result, which agrees with constraints on n from the x-ray cluster 
temperature function, should be roughly independent of the value of O or A. 

Subject headings: Galaxies: clusters: general — cosmology: theory — X-rays: galaxies 



1. INTRODUCTION 

The dynamical state of a system of galaxies may be char- 
acterized by three intimately related physical quantities: 
the bolometric x-ray luminosity Lx, the average emission- 
weighted plasma temperature (T), and the average dark 
matter velocity dispersion, a& m (as traced by the projected 
velocity dispersion of the galaxies, a p ). There is an ob- 
served correlation among these parameters, crudely given 
by (T) oc L x oc (T)" 2 , and L x oc cr° 3 , and roughly 
in agreement with the predictions of physical models. My 
goal in this letter is to show that the value of the slope 
c*3 at z — can constrain the spectrum of the primordial 
density fluctuations, P(k) oc k n . 

One powerful tool for linking the x-ray properties of 
galaxy systems with the cosmological parameters already 
exists. The cluster temperature function, which is directly 
related to the cluster mass distribution, can be used at zero 
redshift to measure n. Henry & Arnaud (1991) analyze 
data from the Einstein satellite to obtain n « —1.7 ±0.55; 
data from the ROSAT and ASCA missions (Markevitch 
1998), corrected for the presence of cooling flows, implies 
a somewhat steeper n « — 2 ±0.3. The redshift evolution of 
the cluster temperature function can break the degeneracy 
between the density parameter, Q, and the normalization 
of the primordial spectrum (Henry 1997). 

Here I introduce a new zero-redshift method. I suggest 
that (X3, the slope of the Lx — <J P relation, is related to 
n in a way that is complementary to the dependence of 
the temperature function on n. While it is important to 
have a flux-limited sample for the temperature function 
method, using 013 to compute n requires only a large num- 
ber of pointed observations of systems of galaxies. In §|| 
I derive the relationship between 013 and n. In I com- 
pare the results with available observations, and discuss 
possible systematic biases. In §|] I summarize. 

2. DERIVATION 

Here I use some of the definitions in the paper by 
Navarro, Frenk, & White (1997; NFW). My results, how- 
ever, are largely independent of their N-body simulations 
of dark matter halos, and apply to density profiles differ- 
ent from the one they introduce. Eke, Navarro, & Frenk 
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(1998) provide a derivation of the Lx — (T) relation which 
is analogous to part of what follows, but does not focus on 
any cosmological use of a 2 ■ 

2.1. Halo Parameters and Definitions 

A spherically symmetric dark matter halo may be char- 
acterized by T200) the radius which encloses 200 times the 
critical density of the universe, p cr it- The characteristic 
mass, M200, is then 
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Here Hq is the Hubble constant, £1 is the density parame- 
ter, and z is the redshift. The halo's characteristic circular 
velocity is V200 = -y/GA^oo/^oo, or , with the substitution 
of equation (|l|) , 

y 20 o = 10H r 2m Z 1/2 . (4) 
Now consider a halo density profile of the form, 

CT 

Pdm{r) = p cnt 5 c p{ ), (5) 

^200 

where c is the halo concentration, S c is the characteristic 
density, and p(y)is a nonnegative, declining function of y. 
Then equation (fil) requires the following relation between 
S c and c, 
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C(y) = / t 2 p(t)dt. 
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In scale- free cosmologies, the characteristic density also de- 
pends on n, the slope of the initial fluctuation spectrum, 
because c should trace the mean density of the universe at 
the time of the halo's formation. NFW use Press-Schechter 
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theory to examine the exact relationship among 5 C , A/200 , 
and n. They find that the characteristic density should 
be proportional to the mean density of the universe at the 
epoch when the nonlinear mass M* is (1) a fixed and (2) a 
very small fraction of the current halo mass. This implies 
that S c should scale with M and n the same way as the 
average background density scales with M* and n: 



Here A(T) is the bolometric emissivity as a function of the 
local temperature T, n e is the electron number density ni 
is the ion number density and the integral is over the entire 
volume of the system. Now if (1) rij = n e , (2) the system is 
spherically symmetric, and (3) the gas density is related to 
the dark matter distribution by p g as(r) — fpdm(r), where 
/ is the gas mass fraction, then 
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The slope of this scaling relation is better than 5% ac- 
curate for all Q < 1 cosmologies with either A = or 
Q + A = 1, where A is the cosmological constant in units 
of 3H 2 . With the substitution 7 = — (n + 3)/2 and use of 
equations (0), (|), and (§), 



— — oc V 3 ~< Z~ l12 
C(c) 200 ' 



(9) 



Next, I compute the circular velocity profile of the 
halo, Vdrc = \/GM(r)/r. Under the transformation 
r -> 2/r 20 o/c, 



V clrc (y) 2 = V, 
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Note that the dark matter velocity dispersion profile, 
(7dm (y), will not in general have the same shape as V C i rc (y), 
contrary to the the assumption in similar derivations (e.g. 
Eke et al. 1998). Rather, the velocity dispersion of the 
halo is given by the Jeans equation for a spherical, nonro- 
tating system of collisionless particles (Binney & Tremaine 
1987): 



dr 



2(3 3 pa 2 r = - 



GMp 



(11) 



Here oy is the radial velocity dispersion, /3j is the velocity 
anisotropy parameter, and M is the mass inside the radius 
r. For systems with isotropic velocity dispersion tensors, 
(3j = 0, and locally a\ m — So 2 . In terms of the previously 
defined quantities, equation (pd|) has the solution 
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Finally, consider the average dark matter velocity disper- 
sion, (<Tj m ), within the virial radius r2oo: 

,2 N Jo 3v 2 (y)p(y)y 2 dy 

( ff dm) = r c -, : 9 , (14) 
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D(y) = / 3a 2 (t)p(t)t 2 dt. 
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2.2. X-Ray Luminosity 
The x-ray luminosity of a ball of plasma is 
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where p = p g as/{nim p ) is the mean molecular weight. 
Substitution of equation dq) yields, 
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Now eliminate r2oo, S c , and p cr u using equations (Q), 
and(|6|): 
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Next, note that thermal bremsstrahlung is the domi- 
nant cooling process for rich clusters of galaxies; hence 
A(T) cx T 1 ! 2 . Then, if the gas is in local hydrostatic equi- 
librium with the dark matter, T(r) oc Ud m (r) 2 , and hence 
A(T, r) cx <7d m (?"). Leaving out the constants and substi- 
tuting equation ([12) yields, 
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The integral on the right-hand side is just a number, and 
may be dropped. 

Now one must replace V200J which is not observable, 
with the projected galaxy velocity dispersion a p , which 
can be determined from optical surveys. In the cores of 
relaxed clusters a 2 is proportional to (<Tj ), the average 



dark matter velocity dispersion (equation 14). Then 
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If / and c remain constant, the above expression reduces to 
the traditional Lx oc a 4 scaling law from simpler, dimen- 
sional arguments (e.g., Quintana & Melnick 1982). How- 
ever, equation (^) tells us that c is a function of the veloc- 
ity dispersion and the slope of the primordial spectrum. 
In fact, once p(y) is specified, equations (R), (||), (p5h, and 
( |l6| ) may be used to eliminate V200 > c, C (c) , and D{c) , and 
one may thus obtain the dependence of Lx on o~ p and n. 
Specifically, suppose that within the range of interest for 
c, C(c) and D(c) exhibit a power law behavior of the form 
C(c) oc c p and D(c) oc c q . Then 
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2.3. Model Dependency 

I consider dark matter density profiles of the form 
p{y) = y a (^ + y b )~ d - Thus each profile's shape may be 
specified by a set of three numbers, (a, b, <£). Some com- 
mon profiles and their properties are listed in Table flj. 

Once the set (a,b,d) is specified, C(c) and D(c) are 
readily computable. The relevant range of c for systems of 
galaxies comes from measurements of surface number den- 
sity profiles and of r2oo in clusters and groups of galaxies 
(Carlberg et al. 1997; Mahdavi et al. 1999). In these works, 
halos with masses in the range 10 14 — 1O 16 M have con- 
centrations c ~ 2.5-9.5, in good agreement with N-body 
simulations (e.g., NFW). For all sets (a,b,d) in Table [l[ 
the power law approximations C(c) cx c p and -D(c) oc c q 
are better than 8% accurate everywhere within c =2.5-9.5. 

Figure [l] shows £ (n) from equation (^4|) for various pro- 
files. In all cases, £,(n) is positive, and approaches zero 
as n — » —3. As n — > 0, the models all predict a signif- 
icant flattening of the Lx — o p relation. This is under- 
standable through equation (||) : the characteristic density 
5 C is highly anticorrelated with halo mass in the n rj 
universes, and hence the emission measure increases quite 
slowly or not at all with mass. As n — » —3, S c becomes 
nearly independent of the halo mass, and therefore Lx 
increases rapidly with the velocity dispersion. 

3. APPLICATION 

For low-redshift clusters of galaxies, Z(z,il) ks 1 with 
high accuracy. The dependence of the gas mass fraction 
/ on (T) and hence a p is poorly determined. Some, using 
standard analysis, claim it should increase slightly with 
(T) (e.g., David, Jones, & Forman 1995); others say that 
accounting for cooling flows should cause it to decrease 
slightly with (T) (e.g. Allen & Fabian 1998). Mohr, Math- 
iesen, & Evrard (1999), in their detailed study of clusters 
with ROSAT pointings, find that within 1 Mpc / is nearly 
independent of (T). Because of the range of contrasting 
findings, I adopt / oc (T) 0±a2 oc cr£ 0±0 - 4 at the 68% 
confidence level. 

There is more agreement among observers regarding the 
empirical value of 0:3: Quintana & Melnick (1982) have 
o?3 = 4.0±0.7 for data from the Einstein satellite; theirs are 
2-10 keV luminosities, which scale similarly to bolometric 
luminosities for rich clusters. More recently, Mulchaey & 
Zabludoff (1998; MZ98) combined ROSAT observations of 
groups and clusters with deep optical spectroscopy to ob- 
tain Q!3 = 4.29 ± 0.37 for bolometric luminosities. I adopt 
the average value, — 4.15 ± 0.4. 

Assuming normally distributed errors, the adopted val- 
ues of / and «3 constrain £ to be < 1.0 at the one-sided 
90% confidence level, or < 1.93 at the one-sided 99% confi- 
dence level. Figure ([j]) shows the £ = 1.0, 1.93 boundaries. 
To accommodate all the models I consider in this scenario, 
n must be < —2.0 at 90% confidence level, and < —1.1 at 
the 99% confidence level. 

At least two systematic effects could bias these results. 
Preheating of the plasma in k(T) 4 keV clusters (Pon- 
man, Cannon, & Navarro 1999) might suppress their lu- 
minosities, while leaving those of k(T) <; 4 keV clusters 
unchanged. One can avoid this bias by ignoring all clus- 
ters with k(T) < 4 keV. I find that removing these clusters, 
which make up « 25% of the MZ98 sample, does not sig- 



nificantly affect the slope of the observed Lx — <J relation. 

Also, cooling flows might affect the luminosities. While 
Markevitch (1998) finds that the slope of the L x - (T) rela- 
tion does not change as a result of removing cooling flows, 
Allen & Fabian (1998) find that a 2 = 3.1 ± 0.6 changes to 
a2 = 2.3 ± 0.4 after including cooling flows. The Marke- 
vitch (1998) method, which removes the cooling flow alto- 
gether, is more appropriate for probing the gravitational 
potential than the Allen & Fabian (1998) method, which 
includes the luminosity of the cooling component. 

To constrain the effect of cooling flows on , I conduct 
the following test. Of the MZ98 clusters, 34% are are also 
contained in Markevitch (1998). I fit 0:3 for just these clus- 
ters, using the Lx which excludes the cooling component. 
I find that the slope changes to 0:3 = 3.85 ± 0.3; with this, 
the upper limits on n become n < —1.7 and n < —0.9 at 
the one-sided 90% and 99% confidence levels. The chief 
result of this paper — that if n were much greater than —1, 
013 should be « 2 instead of the observed value, w 4 — is 
therefore not affected by the 10% correction due to cooling 
flows. 
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Fig. 1. — The correction to the slope of the Ly^ — a v relation, £(n), 
for various density profiles. The dashed lines represent the one-sided 
90% and 99% upper bounds on £ from available observations. The 
arrows indicate the same confidence intervals on n. 

4. CONCLUSION 

If the characteristic densities of clusters of galaxies trace 
the background density of the universe at the time of each 
cluster's formation, there should be a relatively simple re- 
lationship between the x-ray luminosity Lx, the observed 
velocity dispersion tr p , and the slope of the primordial 
power spectrum, n. This relationship, given by equations 
d23| ) and (|^), depends slightly on the density profile of 
the clustered dark matter, but should not significantly de- 
pend on or A when applied to low-redshift clusters of 
galaxies. For a wide range of assumed density profiles, the 
observations imply n < —2.0 and n < —1.1 at the one- 
sided 90% and 99% confidence levels, respectively. This is 
consistent with the bounds from the cluster temperature 



4 Cosmological Lx — cr Relation 

Table 1 
Density Profiles Considered 



Model a,b,d p q 



Hernquist (1990) 1,1,3 0.35 0.16 

NFW 1,1,2 0.74 0.58 

Plummer Sphere 0,2,5/2 0.14 0.05 

King (1962) 0,2,3/2 0.74 0.59 

Jaffe (1983) 2,1,2 0.17 0.03 



function, which give n between —2.3 and —1.15. Improv- 
ing this constraint depends largely on a better understand- 
ing of preheating, cooling flows, and the variation of the 
gas mass fraction with a p . 
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nian Institution and by the National Science Foundation. 
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